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ABSTRACT 

We present phase resolved optical photometry and spectroscopy of the accreting mil- 
lisecond pulsar HETE J1900. 1—2455. Our i?-band light curves exhibit a sinusoidal 
modulation, at close to the orbital period, which we initially attributed to X-ray heat- 
ing of the irradiated face of the secondary star. However, further analysis reveals that 
the source of the modulation is more likely due to superhumps caused by a precessing 
accretion disc. Doppler tomography of a broad Ha emission line reveals an emission 
ring, consistent with that expected from an accretion disc. Using the velocity of the 
emission ring as an estimate for the projected outer disc velocity, we constrain the 
maximum projected velocity of the secondary to be 200 km s~^, placing a lower limit 
of 0.05 Mq on the secondary mass. For a 1.4 Mq primary, this implies that the orbital 
inclination is low, < 20°. Utilizing the observed relationship between the secondary 
mass and orbital period in short period cataclysmic variables, we estimate the sec- 
ondary mass to be '^0.085 Mq, which implies an upper limit of ~2.4 Mq for the 
primary mass. 

Key words: accretion, accretion discs - techniques: photometric - techniques: spec- 
troscopic - pulsars: individual: HETE J1900. 1—2455 - X-rays: binaries 



1 INTRODUCTION 

Low mass X-ray binaries (LMXBs) are close binary systems 
consisting of a degenerate primary - black hole or neutron 
star (NS) - accreting matter from a low mass (< 1 Mq) 
secondary, via Roche lobe overflow. X-ray transients (XRT) 
are LMXBs which undergo peri odic outbursts, explained by 
the d isc instability model (e.g. lOubus. Hameurv fc Lasotal 
I2OOII ). while persistently bright systems are permanently 
in the outbursting state. One subclass of LMXBs are 
the accretion-powered millisecond pulsars (AMSPs), the 
first of which, SAX J1808 . 4-36 58, was discovered in 1998 
IWiinands fc Van der Kii3ll998l ). Seven m ore of the s e sys- 
tems have been discovered since then ("Poutanenl I2OO6I : 
IWiinands 2006; Markwardt . Krimm fc Swank 2007). 

In general, the outburst duration for AMSPs ranges 
from a few weeks to months. AMSPs generally contain a 
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weakly magnetized (~10*-10® G) NS with spin frequencies 
between 180 and 600 Hz. The orbital periods range between 
40 min and 5 h. The second ary star in these sys tems is ei- 
ther a white or brown dwarf llFalanga et al.ll2007l . and refer- 
ences therein). iFalanga et al. I l|2005l ) have measured the pre- 
dicted decrease in the NS spin period for the AMSP system 
ICR J00291-I-5934, supporting the idea that these AMSPs 
are in fact old NS, which have over time been spun up to 
millisecond periods by acquiring angular momentum from 
the accretion of material from the secondary. As such, they 
provide the missing link between LMXBs and old, isolated 
millisecond radio pulsars. 



In comparison to brighter LMXBs, relatively little is 
known about the optical properties of AMSPs. Even in out- 
burst, the eight so far discovered are comparatively faint, 
the brightest bein g SAX J1808.4-3658 at i? ~ 16.2 mag 
IWang et al.ll200ll '). In quiescence, these systems are gener- 
ally very faint optically, the brightest measured again being 
SAX J1808.4-3658 at ~ 20.9 mag. 
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HETE J1900. 1- 2 455 is the seventh AMSP, discovered 
by IVanderspek et alj (|2005l ') on 2005 June 14 by the High 
Energy Transient Explorer II {HETE II). Several type-I 
bursts have since been observed, and assuming that the peak 
flux observed during the brightest of these is Eddington- 
limited, HETE J19 00. 1-2455 is at a distance of ~4.3 kpc 
jSuzuki et al.ll200'i1 ). 

The pulsation frequency was determined with 
the Rossi X-ray Transient Explor er to be 377.3 Hz 
jMorgan. Kaaret fc Vanderspekl |2005| ) . Pulse timing anal- 
ysis by iKaaret et al. I 12006') revealed a circular orbit with 
a period (Porb) of ~83.3 minutes (4995.258 ± 0.005 s), 
and a projected primary semimajor axis, aisini, of 18.41 
± 0.01 X 10~^ light second. These parameters combine to 
provide a mass f unctio n /(M) = 2.004 ± 0.003 x 10"** Mq. 
iGallowav et al.l (|2006l ) report that pulses were detected 
intermittently in the first two months after discovery. No 
pulses have been reported since then. 

Unusually for an XRT, this system has remained active 
for more than 2 years, at a level of ~2 x 10"^® erg s~^. Several 
times during 2007 March/February, the X-ray flux fell by an 
order of magnitude, but each time return ed to outburst lev- 
els within ~1 week (|Gallowav et al.ll2007l '). In 2007 May, the 
source was observed to decline to its lowest level since discov- 
ery (by a factor of > 10"^), lasting for several weeks, prompt- 
ing speculatio n that HETE J1900. 1—2455 was returnin; 
to quiescenc e (Pegena ar et al. I l2007al : [GaTlowav et al.|[200" 
[Torres et al.l [2007). However, the system returned to its 
outburst state within ~2-3 weeks l|Degenaar et al ] |2007bl : 
iGarnavich et al.ll2007l ). 

The optical counterpart was found by|^ (|2005l ) with 
an P-ba nd magnitude of ~ 18.4 mag. Subsequent observa- 
tions bv lSteeghs et aL I l|2005ll found the i?-band magnitude 
to be 18.02 ± 0.03 mag, and the V — R colour (dereddened) 
to be —0.16 mag, with spectroscopy revealing a broad He ll 
A4686 A emission line. 

In this paper we present optical photometry and spec- 
troscopy of HETE J1900. 1-2455 in its outburst state. These 
observations were motivated by several goals, key amongst 
them the aim of measuring the radial velocity of the sec- 
ondary, and hence constraining the mass of the compact 
object. For the longer period transient systems, such mea- 
surements are possible via radial velocity studies of ab- 
sorption lines from the secondary once the system has re- 
turned to quiescence. For the shorter period AMSPs how- 
ever, the secondary is likely to be so faint as to make such 
observations extremely di fficult, if not impossible. Ho wever, 
ISteeghs fc CasaresI l|2002l ) and ICasares et ahl l|2004l ') have 
shown how radial velocity measurements of the secondary 
are still possible for X-ray bright systems, by studying the 
Bowen blend emission from the irradiated face of the sec- 
ondary: this technique may offer the only way of measuring 
the velocity of the secondary in AMSPs (i.e. whilst they are 
still X-ray bright), and usefully constraining the mass of the 
neutron star. 



2 DATA 

Our data consist of photometry obtained with the 2.3 m Ad- 
vanced Technology Telescope (ATT) at Siding Spring Obser- 




« 



• 



I 



Figure 1. LBT image of the field of HETE J1900.1-2455 (Sloan 
r-band), 2. '3 x 1.'6. The black arrow indicates the position of 
HETE J1900. 1-2455. 



vatory, Australia, the 3.5 m Wisconsin-Indiana- Yale-NOAO 
(WIYN) telescope at Kitt Peak, Arizona, and the 8.4 m 
Large Binocular Telescope (LBT) at Mount Graham, Ari- 
zona. Our spectroscopy was obtained using the High Res- 
olution Echelle Spectrometer (HIRES) on the 10 m Keck I 
telescope at Mauna Kea, Hawaii. 



2.1 Photometry 

We obtained photometry on 2 nights in 2006 September with 
the 2.3 m ATT, using the CCD4240 e2v detector, mounted 
at the //18 Nasmyth focus. The detector has 2148 x 2148 
13.5 ^.m pixels, with 0."34 pixePV On 2006 September 21 
UT, we obtained 16 _R-band exposures, covering slightly more 
than 1 orbital period, with exposure times of 60 s (x 1), 300 
s (x 9) and 240 s (x 6). On the following night, we obtained 
R, V and _B-band exposures, 2 x 300 s exposures of each. 
On both nights, expos ures were take n of the photometric 
standard field SA-112 (|Landoltj[l99i ). Seeing was ~2" on 
both nights, but deteriorated towards the end of the second 
night. 

Three 300 s i?-band exposures were also taken with the 
3.5 m WIYN telescope on 2006 September 19 ut, using the 
WIYN Mini-Mosaic Imager, consisting of 2 SITe CCDs, each 
with 4096 X 2048 15 ^.m pixels, with 0."14 pixeP^ The see- 
ing was ~2" for all 3 exposures. Expo sures were also taken 
of the standard star field PG031-I-051 (|Landohlll993 '). 

After the decline/rebrightening episode of 2007 
May/June, we obtained additional observations using the 
LBT. On 2007 June 13 ut and June 24 UT we obtained 200 
s Sloan r-band exposures of HETE J1900. 1-2455 with the 
Large Binocular Camera (Bl ue) (Giallongo et al., in prepa- 
ration; (RagazzonTeF^aLjlloO^) at the prime focus of the LBT 
(20 exposures per night). The eev-blue detector consists of 4 
CCDs with 2048 x 4608 13.5 ^m pixels, with 0."224 pbceP^ 
The seeing on 2007 June 13 was ~0. "85-1. "15, while on June 
24 was ~1"-1."5. The HETE J1900.1-2455 field is shown 
in Fig. [U (from LBT on 2007 June 13, with seeing ~0."85). 

All the photometry data were reduced in the same way. 
The frames were firstly bias corrected, trimmed and fiat- 



Optical observations of HETE J1900. 1-2455 3 



fielded using the CCDPROC routines i n IRAI^I. Phot ometry 
was performed using the daophot (|Stetsonl Il992l 'l point 
spread function fitting package in iraf. The magnitude of 
HETE J1900. 1—2455 was measured relative to several bright 
stars. These relative magnitudes were then calibrated by 
comparison with the standard star frames. Note that the 
WIYN data was used for photometric calibration only. 

The ATT data from 2006 September gives R = 17.91 ± 
0.05 mag, V = 18.00 ± 0.04 mag and B = 17.92 ± 0.04 
mag. The errors quoted are an estimate of the systematic 
uncertainty. We note that the me asured fl-band magn itude 
is consistent with that reported bv lSteeghs ilaDiiooi). On 
2007 June 13, R = 18.51 ±0.05 mag, while on 2007 June 24, 
R = 18.44 ± 0.05 mag. 

Taking the value for the hydrogen column den- 
sity to the source, A^h, as 1.6 ± 0.4 x 10^^ cm~'^ 



l|Campana. Cucchiara fc Burrow^ 20051). a nd using the 
method detailed bv lSavage fc MathisI (|l979l ) yields E{B - 
V) = 0.30 ± 0.08. The resulting dereddened colours are 
{B - V)a = -0.34 ± 0.08 mag, and (V - R)a = -0.17 
± 0.06 m ag (again, consistent with the value reported by 
ISteeghs et al. 20051. 

In the i?-band, HETE J1900. 1-2455 is blended with a 
faint star 2" to the northwest (see Fig. [TJ. This star was 
fitted simultaneously with HETE J1900. 1-2455, and was 
found to have an _R-band magnitude of 20.1 ± 0.1 mag, 
which was constant, within these errors, in all our observa- 
tions. This star was not visible in the B or V -band frames. 
This is consistent with the observations of ISteeghs et al.l 
l|2005h . who note that this blended star is ~2 magnitudes 
fainter (in the ii-band) than HETE J1900. 1-2455. 

In the LET data, in frames with good seeing, an ad- 
ditional faint star is observed ~l" northwest of HETE 
J1900. 1—2455, with an _R-band magnitude ~5 mag fainter 
than HETE J1900. 1-2455. 

In Fig. [2] we plot the three _ R-band light curves, folded 
on the pulse timing ephemeris of iKaaret et al.l (|2006l ). The 
phase error at the epochs of these optical observations, ex- 
trapolated from the the error in the period of the pulse tim- 
ing ephemeris, is ~0.01. Each light curve exhibits a similar 
modulation, but that of 2007 June 24 has an amplitude ~3 
times greater than for the other two. Clearly, the phasing 
of the three light curves is not consistent, implying that the 
period of this modulation is different to the orbital period. 



2.2 Spectroscopy 

Our spectroscopy wa s obtained using the HIRES spectrome- 
ter (|Vogt et al.lll994l ') at the right Nasmyth focus of the 10 m 
Keck I telescope, with a 52.68 lines mm~^ echelle grism and 
the HIRESr 250 lines mm~^ long wavelength cross disperser. 
The slit width was 0."86 with a slit length of 7". The echelle 
grating angle was 0° and the cross disperser grating angle 
was 0.°2405, providing a wavelength coverage of ~ AA3900- 
8500 A, spread over a mosaic of 3 MIT-Lincoln Labs CCDs 
("red", "green" and "blue"), each with 2048 x 4096 15 (im 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



00 



(a) 



Pi J I 



1^ - - 1 i 
^ ■ 1 



1 1 



I I 



H 1 1 1 1 \ 1 h 



H 1 1 1 h 



00 



(b) 



lO - I ■ 

o6 - 



CD 

o6 



ii i 



H 1 1 1 1 1 \ 1 1 1 \ 1 1 1 1 1 1- 



o6 



Ln 
od 



CD 

od 



0.5 



1 .5 



Orbital phase 



Figure 2. R-band light cur ves, folded on the pulse timing 
ephemeris l lKaaret et al.|[2006h . and plotted twice for clarity, (a) 
ATT light curve from 2006 September 21, (b) LET light curve 
from 2007 June 13, (c) LET light curve from 2007 June 24. 



pixels, with 0."12 pixel"^ in the spatial dimension. The im- 
ages were binned by 2 pixels in the spatial dimension. 

15 echelle spectra were obtained with HIRES on 2006 
October 11 UT. Appropriate bias, fiat, trace star, ThAr ar- 
clamp and fiux standard (G191B2B) exposures were also 
taken. The object exposure times were 300 s (x 5), 360 s 
(x 8) and 420 s (x 2), with the 15 spectra covering slightly 
more than one orbital period. 

The spectra were extracted using the makeeQ (MAuna 
Kea Echelle Extraction) HIRES reduction package, which 
takes the two dimensional raw fits images, and produces 
optimally extracted, wavelength calibrated, one dimensional 
spectra. The lines of interest to us were on the "blue" and 
"green" CCDs. The wavelength calibration was performed 
using the ThAr arc lamp exposures, by fitting a sixth or- 
der polynomial to an average of 30 lines per echelle order 
(minimum 20); this resulted in a dispersion of 0.015-0.025 
A pixel"^ for the "blue" CCD and 0.021-0.032 A pixel"^ for 
the "green" CCD, and an RMS scatter of < 0.003 A (average 
~0.002 A). To confirm the extraction method implemented 
in MAKEE, the spectra were also extracted using the hires 
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Figure 3. Average spectrum for Ha, normalized to the local 
continuum. 



REDUX softwar^l, and the echelle package in iraf: how- 
ever the spectra extracted using all three methods were very 
similar. 

The main emission features present in our spectra are 
Ha A6563 A, He ll A4686 A and the Bowen blend of N ill 
and C HI lines centred at ~A4640 A. H/3 is present as a 
weak emission line within a broader absorption feature. The 
echelle orders containing these features were exported for 
further processing with the molly spectroscopic data anal- 
ysis package, and doppler, a package used for performing 
maximum entropy Doppler tomography. The spectra were 
re- binned onto a velocity scale of 40 km s~^ pixel" ^, yielding 
a signal-to- noise (S/N) ratio of 6-16 pixel" ^. The S/N varied 
across each order, particularly on the "blue" CCD, with the 
lowest S /N near the edges of the orders. For each spectrum, 
the continuum was fitted with a 2nd order polynomial, af- 
ter masking the emission regions. The fit was then divided 
into the spectrum, and the remaining continuum subtracted. 
The average processed spectrum for Ha is shown in Fig. [S] 
and that for the Bowen blend and He ll A4686 A is in Fig. 
m Note that the data in Fig. 2] has been boxcar smoothed, 
and shifted to the rest frame of the system. From the av- 
eraged spectra, the equivalent width for Ha is —3.1 ± 0.3 
A, for He 11 A4686 A is -2.5 ± 0.3 A, and for the Bowen 
blend is -2.1 ± 0.2 A. This is in contrast to the observa- 
tions of lSteeghs et al.l l|2005l l , who report that although their 
spectroscopy was contaminated by a nearby star, the main 
feature observed was a broad emission line from He ll A4686 
A and Balmer emission was at most very weak. 



3 RESULTS 

3.1 i?-band Light Curves 

The three light curves in Fig. [5] exhibit a sinusoidal modula- 
tion, classically attributed to X-ray irradiation of the inner 
face of the secondary, where the modulation is due to the 
phase dependent visibility of this higher temperature region. 
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Figure 4. Average spectrum for the Bowen blend and He n 
A4686 A, normalized to the local continuum. The spectrum has 
been shifted to the rest frame of the system, based on the sys- 
temic velocity determined from the Ho emission line in section 
3.2. The spectrum has also been boxcar smoothed. 



The phasing of the ATT light curve is approximately con- 
sistent with this scenario, with the brightest measurements 
occurring near phase 0.5. However, the phasing of the other 
light curves is completely different, suggesting that the ob- 
served modulation does not vary on the orbital period. In 
particular, the two LBT light curves are separated by only 
11 nights, and in that time the phasing of the light curves, 
folded on the orbital period, differs by ~0.2. 



3.2 Ha emission 

Fig. [3] shows the average Ha spectrum. This broad emission 
line exhibits a hint of a double peak (a classic signature of an 
accretion disc), although the peaks are not as distinctive as 
those seen in the spectra of other LMXBs. The distribution 
of the flux appears to be symmetrical about a central wave- 
length, suggesting that, on average, the emission from the 
disc is axisymmetrical. To constrain the systemic velocity, 7, 
we fitted a Gaussian to the wings of the averaged spectrum, 
where the emission is from disc regions closer to the cen- 
tre, and presumably less susceptible to corruption fr om the 
hotsp ot and other non-uniformities in the outer disc (iMarshI 
119981 ). This gives a value for 7 of 132 ± 8 km s"^ 

We investigated the emission line distribution of the 
system by examining a Doppler tomogram of the Ha line 
(Fig. [5]), constructed using 7 = 132 km s~^. As with the 
photometry, the spectra were phased using the pulse timing 
ephemeris. For this work, we used the maximum entropy 
method, implemented in do ppler: see iMaJsh and Horng 
(1988), Marsh (200i|) andHteeghs (200i) for further details 
and applications of Doppler tomography. 

Enhanced emission is present in the { — V^,+Vy) quad- 
rant, extending into the (-V^,-Vy) quadrant, with a faint 
ring of emission. Such a ring of emission is seen in Doppler 
tomograms of many LMXBs, and is attributed to emission 
from the outer disc (iMarsbl [l998l ). This emission ring lies 
between radii of ~200-400 km s~^. Plotting only the ax- 
isymmetric component of this tomogram reveals that the 
brightest emission is from an annulus at a velocity of ~230 



Optical observations of HETE J1900. 1-2455 5 




Figure 5. Doppler tomogram for Ha. The Roche lobe of the sec- 
ondary is plotted. The lower curve represents the gas stream ve- 
locity (the free-fall ballistic trajectory from the inner Lagrangian 
point, Li) and the upper curve is the Keplerian velocity of the 
accretion disc along the gas stream. Both curves are marked with 
a circle when the gas stream moves O.IRli (0.1 of the distance be- 
tween the primary and Li) closer to the primary. The gas stream 
trajectory and disc velocity along the stream were calculated us- 
ing K2 = 200 km s-i and q = 0.035. 



km s^^, but with significant emission extending to ~150 
km s~^. Since the secondary must be moving slower than 
the outer disc, we choose 200 km s~^ as an upper limit to 
K2, the velocity of the secondary projected onto the line of 
sight. 

For illustration, we plot in Fig. [S]the Roche lobe of the 
secondary, the gas stream trajectory and accretion disc ve- 
locity curve, computed using K2 = 200 km s~^. Since Ki 
is already accurately known from the pulse timing analy- 
sis (2 7raisint/Porb ~ 6.942 ± 0.003 km s~^; iKaaret et al.l 
I2OO6I '). q (= Ki/K2)= 0.035. K2 and q are the only param- 
eters upon which the locations of the Roche lobe and the 
gas stream trajectory depend. At this value of K2, the pre- 
dicted gas stream misses the hotspot location. Lower val- 
ues of K2 would cause the gas stream to be closer to the 
hotspot. However, we also note that Doppler tomograms of 
the He 11 A4686 A emission line in the two NS LMXBs V801 
Ara and V926 Sco also show the peak of emission in the 
{ — V^, — Vy) quadrant, close to the —V^ axis, and for these 
systems this emission is attributed to an extended disc bulge 
l|Casares et al.|[2006l ). As will be discussed in section 4, the 
location of the hotspot in HETE J1900. 1-2455 may be ex- 
plained by a warped, irradiated accretion disc. 



3.2.1 Projected secondary velocity vs. inclination 

Using the standard mass function equation, and with K2 = 
Ki/q — 27raisim/gPorb it is possible to calculate q and K2 
for various orbital inclinations (i) and primary masses (Mi). 




i n 

Figure 6. Projected secondary velocity {K2) vs. orbital inclina- 
tion angle (i) for a primary mass of 1.4 Mq. The upper limit on 
K2 is set by the velocity of the Ho emission ring. 



Fig.|6]shows a plot of K2 vs. i for a primary mass of 1.4 Mq. 
If K2 < 200 km s-\ i < 20°, for all Mi > 1.4 Mq. 

Hence our upper limit on K2 provides an upper limit 
on i of ~20° . This low inclination is consistent with the fact 
that the double peak in Ha is not as clear as in systems with 
higher inclinations, and thus higher projected disc velocities 
(see for example the LMXB XTE J1118-f480, with i ~70°; 
lElebert. Callanan fc Torres, 2OO61 '). 



3.3 Bowen blend emission 

ISteeghs fc CasaresI (|2002l ) discovered that emission from the 
Bowen blend (a blend of N lll, C lll and O 11 lines, centred 
at ~A4640 A) traced the motion of the irradiated secondary 
star in the LMXB Sco X-1. Some of the major features 
seen in the Bowen blend in X-ray binaries are the N ill 
lines at A4634.13 A, A4640.64 A and AA4641. 85/4641. 96 
A, and the C ill lines at A4647.4 A, A4650.1 A 
and A4651 A (j Hvnes et al' '200l"; ISteeghs fc CasaresI I2OO2I : 
iNelemans. Jonk cr & Stccghs 2006). Our average spectrum 
of the Bowen blend reveals a broad range of emission be- 
tween A4620 A and A4650 A (Fig. |3| . Superimposed on this 
broad emission, there appear to be several narrower emission 
features, three of which correspond to the expected wave- 
lengths of the N III A4634.13 A line, the N ill blend near 
A4641 A, and the A4647.4 A C ill line. 

Because of the low S/N, we performed Doppler tomog- 
raphy of all of the lines in the Bowen blend simultaneously, 
with the intention of finding emission from the irradiated 
face of the secondary. Although this map contains a sin- 
gle sharp spot of emission, it is in the {+Vx,+Vy) quadrant, 
~45° away from the expected location of secondary emission 
on the +Vy axis. The significance of this emission is uncer- 
tain because of the low S/N of the individual exposures. 
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3.4 He II A4686 A emission 

The He ll A4686 A emission is contaminated in the wings, 
particularly on the red side. While the lower S/N in the 
wings might explain some of this variability, is does not ac- 
count for all of the systematic artefacts seen - in particular, 
the 2 dips on the red side, near A4704 A and A4710 A, can 
be seen in averages of just a few spectra, and so are not due 
to random noise. As well as this, the average line appears 
asymmetrical, with an excess of flux between AA4695-4705 
A. The S/N ratio throughout the order is also lower than 
that for Hq. Despite this, a Gaussian fit to the wings of this 
line (masking the line core, and the excess on the red side) 
gives a systemic velocity consistent with that seen in the Ha 
emission. 



4 DISCUSSION 

With its maximum brightness near phase 0.5, we attempted 
to fit the first of our light curves (Fig.[2fa)) using the Eclips- 
ing Light Curve (elc) code (|Orosz fc Hauschildt.2000, '). un- 
der the assumption that the modulation was due to X-ray 
heating. This provided a satisfactory fit, although because 
of the degeneracy between inclination angle, disc tempera- 
ture and disc flaring angle, it was not possible to extract a 
unique parameter set, and a wide range of parameter values 
yielded equally good flts. 

However, the phasing of the subsequent light curves 
(Fig. [2Jb),(c)) suggests that the observed modulation is not 
due to X-ray heating of the secondary: the most likely ex- 
planation for this is that the observed modulation is due to 
superhumps (see below). Therefore the apparent alignment 
of the first light curve is most likely coincidental, and the 
ELC code (which does not incorporate superhumps) cannot 
be used to extract reliable system parameters. 

Superhumps are periodic optical modulations, origi- 
nally seen during superoutbursts of SU UMa dwarf no- 
vae (DNe), but also seen in some LMXBs in outburst 
l|Haswell et al.l l2001al . and references therein). The phe- 
nomenon has also been seen in persistently bright systems. 
The superhump period is typically a few percent longer 
than Porb- The current understanding is that in systems 
where the mass ratio is less than ^^0.25, the tidal trunca- 
tion radius lies outside the 3:1 resonance radius. When the 
accretion disc radius extends out to beyond the 3:1 reso- 
nance radius, the disc becomes eccentric, and the perturba- 
tion from the secondary ca uses a prograde prece ssion. For 
precessing discs in LMXBs, iHaswell et al.l (|2001al ) proposed 
that the disc area changes on the superhump period, caus- 
ing the area of the disc intercepting the X-rays to increase. 
The increased reprocessing causes the disc to brighten, giv- 
ing rise to the observed superhumps in LMXBs. Several 
systems exhibit superhumps with almost exact sinusoidal 
modulations like those seen in HETE J1900 .1— 2455, in par- 
ticula r another LMXB, XTE J1118-I-480 (jPatterson et all 
I2OO5I ). The peak-to-peak amplitude of the modulation seen 
in HETE J1900. 1-2455 was ~0.1 mag on 2006 September 
21 and 2007 June 13, but increased to ~0.3 mag on 2007 
June 2 4. The latter amplitud e is quite high, but as pointed 
out bv lHaswell et al.1 (|2001bl ) superhumps due to irradiated 
discs are likely to appear more prominent at low inclinations: 



hence the large superhump amplitude supports our estimate 
of the inclination in section 3.2.1. For an X-ray luminosity 
of ~10''^ erg s~^ the di sc is likely to also su ffer irradiation- 
driven warping (Foulkes. Haswell fc Murravi.2006 . in partic- 
ular their model number 1) and the disc shape continuously 
changes under the combination of prograde apsidal preces- 
sion and retrograde warped precession. Such models may 
also explain the factor of ~3 increase in modulation am- 
plitude over a relatively short timescale (Haswell, private 
communication) . 

The location of the hotspot in the Ha Doppler tomo- 
gram may be interpreted in terms of this irradiated pre- 
cessing disc model: as the warped disc precesses, the gas 
stream/accretion disc i mpact point varies. How this occurs 
can be seen clearly in the lFoulkes. Haswell fc Murravl (|2006| ) 
simulation^ These simulations suggest that multi-epoch 
Doppler tomography of HETE J1900. 1-2455 should reveal 
a hotspot whose location varies between the ( — Vx,+Vy) and 
(-14,-1^) quadrants. 

Previous work has shown that, once the superhump 
period (-Psh) is identified, then the period excess, e (= 
(Psh — Porb) /Porb), cau be dlrcctly related to the mass ratio 
( Patterson et al.ll2005l ). A Fourier transform of our photom- 
etry, however, does not allow us to identify a unique super- 
hump period, because our sampling is too sparse. A more 
comprehensive photometry campaign is required to measure 
Psh, and estimate q. 

Using the empirical relationship between the abso- 
lute visual magnitude of an irradiated accretion disc 
(Mv), the X-ray lu minosity of the central source 
and orbital period (Ivan P aradiis fc M cClintoc^ I1994I : 
Ide Jong, van Paradiis fc Augusteiin 499(^). we find My = 
3.6 ± 0.5 mag. At a distance of 4.3 kpc (jSuzuki et al.ll2007l '). 
such a disc should have an apparent magnitude of 16.7 ± 
0.5 mag. This is consistent with the measured l/-band mag- 
nitude of ~17.1 mag, implying that the optical light, in the 
l/-band, is dominated by an X-ray heated accretion disc. 
Thus, X-ray heating of the disc can explain the observed 
optical b rightness of HE TE J1900.1-2455. 

Kaaret et all (|2006l ') estimated the mass of the sec- 
ondary to be between 0.016 and 0.07 Mq, for primary 
masses between 1.4 and 2.2 Mq, but this was based on an 
assumption of inclination angle: 90° for the lower mass limit, 
and a uniform a priori distribution for the 95 per cent confi- 
dence upper limit. We have seen that for K2 < 200 km s~^, 
i < 20° (see Fig. E)). At this upper limit oi K2, q = 0.035, 
and the minimum value of M2 (for Mi = 1.4 Mq) is there- 
fore 0;05M£7i^_^ 

[Patterson et al.l (|2003l ) have computed a mass-radius 
relation for the secondaries in short period cataclysmic 
variables (CVs), finding that in general, the mass-radius 
curve lies above that for zero-age main sequence stars. 
An empirical equation relates M2 to Porb: M2 = 
0.0764Po;b''a~^°'^ Mq, where a is the factor relating the 
radius of the secondaries in these CVs to that of main se- 
quence stars w i th the same mass, and Porb is in hours. 
[Patterson et all (|2003l ') find that a ~ 1.18. Assuming that 
the secondary in HETE J1900. 1 — 2455 is similar to the sec- 
ondary in these CV systems, this equation predicts M2 = 
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0.085 Mq. With the lower hmit on q of 0.035, this sets an 
upper hmit to the primary mass of Mi < 2.4 Mq. The sec- 
ondary may also be a brown dwarf, in which case Ad2 is 
poorly constrained, as is A/i. 

Finally, we briefly examine the radial velocity of 
HETE J1900. 1-2455 relative to the local Galactic radial 
velocity. Based on its position and estimated distance, 
HETE J1900. 1-2455 is ~3.6 kpc from the Galactic centre, 
~1 kpc belo w the Ga l actic plane. Using the Galactic rota- 
tion curve of lClemen^ (Il985l ). we find that the radial velocity 
of HETE J1900. 1-2455 is ~90 km s"^ relative to the lo- 
cal Galactic rotation. Combined with future proper motion 
studies, this velocity may be used to derive the space veloc- 
ity of HETE J1900. 1-2 455 and compute its orbit around 
the Galactic Centre (e.g. iMirabel fc Rodriguesll2003l) . 



5 CONCLUSIONS 

We have presented a detailed optical study of the AMSP 
HETE J1900. 1-2455. The ij-band light curves reveal a 
modulation whose phase is inconsistent with simple heat- 
ing of the secondary star /accretion disc: we interpret this as 
due to superhumps associated with a precessing accretion 
disc. Fitting to the wings of the averaged Ha emission line 
gives the systemic velocity of 132 ± 8 km s~^, and a similar 
value is found from the He ll line. 

Doppler tomography of the Hq emission line suggests 
a upper limit of ~200 km s~^ on the projected secondary 
velocity, implying that HETE J1900. 1-2455 has a low incli- 
nation, < 20°, and a minimum secondary mass of 0.05 Mq. 

Assuming that the secondary in HETE J1900. 1-2455 
is similar to those seen in short period CVs, we estimate a 
secondary mass of 0.085 Mq, which implies an upper limit 
to the NS mass of 2.4 Mq (again for K2 < 200 km s~^). 
If the secondary is a brown dwarf, then this constraint is 
relaxed. 



ACKNOWLEDGMENTS 

We acknowledge the use of the 2.3 m Advanced Technology 
Telescope of the Australian National University at Siding 
Spring Observatory. The WIYN Observatory is a joint fa- 
cility of the University of Wisconsin-Madison, Indiana Uni- 
versity, Yale University, and the National Optical Astron- 
omy Observatory. Some of the data presented herein were 
acquired using the Large Binocular Telescope (LBT). The 
LBT is an international collaboration among institutions in 
the United States, Italy and Germany. LBT Corporation 
partners are: The University of Arizona on behalf of the 
Arizona university system; Istituto Nazionale di Astrofisica, 
Italy; LBT Beteiligungsgesellschaft, Germany, representing 
the Max-Planck Society, the Astrophysical Institute Pots- 
dam, and Heidelberg University; The Ohio State University, 
and The Research Corporation, on behalf of The Univer- 
sity of Notre Dame, University of Minnesota and University 
of Virginia. We thank LBT director R. Green for alloca- 
tion of DD time, and D. Thompson, LBTO for acquiring 
some of the data. The spectroscopic data presented herein 
were obtained at the W. M. Keck Observatory, which is op- 
erated as a scientific partnership among the California In- 



stitute of Technology, the University of California and the 
National Aeronautics and Space Administration. The ob- 
servatory was made possible by the generous financial sup- 
port of the W. M. Keck Foundation. This research made 
use of NASA's Astrophysics Data System, and the SIM- 
BAD database, operated at CDS, Strasbourg, France. X- 
ray quick-look results provided by the ASM/RXTE team. 
We acknowledge the use of molly and doppler software 
packages developed by T. R. Marsh, University of Warwick. 
We thank J. X. Prochaska, UCO/Lick Observatory, for re- 
ducing our Keck spectroscopy using hires redux. PE and 
PJC acknowledge support from Science Foundation Ireland. 



REFERENCES 

Campana S., Cucchiara A., Burrows D. N., 2005, ATel, 543 
Casares J., Steeghs D., Hynes R. I., Charles P. A., Cor- 
nelisse R., O'Brien K., 2004, RevMexAA Conference Se- 
ries, 20, 21 

Casares J., Cornelisse R., Steeghs D., Charles P. A., Hynes 
R. I., O'Brien, K., Strohmayer, T. E., 2006, MNRAS, 373, 
1235 

Clemens, D. P, 1985, ApJ, 295, 422 

Degenaar N., Wijnands R., Galloway D. K., Lewin 
W. H. G., Homan J., Chakrabarty D., Campana S., Miller 
J. M., Cackett E. M., 2007a, ATel, 1091 

Degenaar N., Campana S., Galloway D. K., Lewin 
W. H. G., Homan J., Chakrabarty D., Morgan E. H., 
Jonker P. G., Cackett E. M., Miller J. M., Wijnands R., 
2007b, ATel, 1106 

de Jong J. A., van Paradijs J., Augusteijn T., 1996, A&A, 
314, 484 

Dubus G., Hameury J.-M., Lasota J.-P., 2001, A&A, 373, 
251 

Elebert P., Callanan P. J., Torres M. A. P., 2006, in Meurs 
E. J. A., Fabbiano G., Eds., Proceedings of the 230th Sym- 
posium of the International Astronomical Union, Cam- 
bridge University Press, pp. 57-58 

Falanga M., Kuiper L., Poutanen J., Bonning E. W., Di 
Salvo T., Goldoni P., Goldwurm A., Shaw S. E., Stella L., 

2005, A&A, 444, 15 

Falanga M., Poutanen J., Bonning E. W., Kuiper L., 
Bonnet-Bidaud J. M., Goldwurm A., Hermsen W., Stella 
L., 2007, A&A, 464, 1069 

Foulkes S. B., Haswell C. A., Murray J. R., 2006, MNRAS, 
366, 1399 

Fox D. B., 2005, ATel, 526 

Galloway D. K., Morgan E. H., Krauss M. I., Kaaret P., 

Chakrabarty D., 2006, ApJ, 654, L73 
Galloway D. K., Morgan E. H., Chakrabarty D., Kaaret P., 

2007, ATel, 1086 
Garnavich P. M., Callanan P. J., Elebert P., Reynolds 

M. T., 2007, ATel, 1110 
Haswell C. A., King A. R., Murray J. R., Charles P. A., 

2001a, MNRAS, 321, 475 
Haswell C. A., Rolfe D. J., King A. R., Murray J. R., 

Charles P. A., 2001b, ApSSS, 276, 41 
Hynes R. I., Charles P. A., Haswell C. A., Casares J., Zurita 

C, Serra-Ricart M., 2001, MNRAS, 324, 180 
Kaaret P., Morgan E. H., Vanderspek R., Tomsick J. A., 

2006, ApJ, 639, 963 



8 P. EleheH et al. 



Landolt A. U., 1992, ApJ, 104, 340 

Markwardt C. B., Krimm H. A., Swank J. H., 2007, ATel, 
1108 

Marsh T. R., Home K., 1988, MNRAS, 236, 269 

Marsh T. R., 1998, in HowcU S., Kuulkers E., Woodward 

C, Eds., ASP Conf. Ser. Vol. 137, Wild Stars in the Old 

West, ASP, San Francisco, pp. 236-244 
Marsh T. R., 2001, in Boffin H. M. J., Steeghs D., Cuypers 

J., Eds., Astrotomography: Indirect Imaging Methods 

in Observational Astronomy, LNP Series 573, Springer, 

Berlin, pp. 1-26 
Miirabel I. F., Rodrigues I., 2003, A&A, 398, L25 
Morgan E. H., Kaarot P., Vandcrspok R., 2005, ATcI, 523 
Nelemans G, Jonkcr P. G., Stccghs, D., 2006, MNRAS, 370, 

255 

Orosz J. A., Hauschildt P. H., 2000, A&A, 364, 265 

Patterson J., Thorstensen J. R., Kemp J., Skillman D. R., 
Vanmunster T., Harvey D. A., Fried R. A., Jensen 
L., Cook L. M., Rea R., Monard B., McCormick J., 
Vclthuis F., Walker S., Martin B., Bolt G., Pavlcnko E., 
O'Donoghue D., Gunn J., Novak R., Masi G., Garradd 
G., Butterworth N., Krajci T., Foote J., Beshore E., 2003, 
PASP, 115, 1308 

Patterson J., Kemp J., Harvey D. A., Fried R. E., Rea R., 
Monard B., Cook L. M., Skillman D. R., Vanmunster T., 
Bolt G., Armstrong E., McCormick J., Krajci T., Jensen 
L., Gunn J., Butterworth N., Foote J., Bos M., Masi G., 
Warhurst P., 2005, PASP, 117, 1204 

Poutanen J., 2006, AdSpR, 38, 2697 

Ragazzoni R., Giallongo E., Pasian F., Baruffolo A., 

Bertram R., Diolaiti E., Di Paola A., Farinato J., Gen- 
tile G., Hill J., Lombini M., Pcdichini F., Spcziali R., 
Smarcgha R., Vernet E., 2006, in Stcpp, L. M., Ed., Soci- 
ety of Photo-Optical Instrumentation Engineers (SPIE) 
Conference, Volume 6267, Ground-based and Airborne 
Telescopes 

Savage B. D., Mathis J. S., 1979, ARA&A, 17, 73 

Steeghs D., 2004, AN, 325, 185 

Steeghs D., Casares J., 2002, ApJ, 568, 273 

Steeghs D., Torres M. A. P., Garcia M. R., McClintock 
J. E., Miller J. M., Jonker P. G., Callanan P. J., Zhao P., 
Berlind P., Hutchins R., Watson C, 2005, ATel, 543 

Stetson P. B., 1992, in Worrall D. M., Biemesderfer C, 
Barnes J., Eds., Astronomical Data Analysis Software and 
Systems I, ASP Conf. Ser., 25, 297, pp. 297-306 

Suzuki M., Kawai N., Tamagawa T., Yoshida A., Naka- 
gawa Y. E., Tanaka K., Shirasaki Y., Matsuoka M., Ricker 

G. R., Vanderspek R., Butler N., Lamb D. Q., Graziani 

C, Pizzichini G., Sato R., Arimoto M., Kotoku J., Maetou 
M., Yaniauchi M., 2007, PAS J, 59, 263 

Torres M. A. P., Rodriguez-Gil P., Steeghs D., Corral- 
Santana J. M., Casares J., Jonker P. G., 2007, ATel, 1090 

Vanderspek R., Morgan E. H., Crew G., Graziani C, 
Suzuki M., 2005, ATel, 516 

van Paradijs J., McChntock J. E., 1994, A&A, 290, 133 

Vogt S. S., Allen S. L., Bigelow B. C, Brosoo L., Brown B., 
Cantrall T., Conrad A., Couture M., Dclancy C, Epps 

H. W., Hilyard D., Hilyard D. F., Horn E., Jcrn N., Kanto 

D. , Koane M. J., Kibrick R. L, Lewis J. W., Osborne J., 
Pardeilhan G. H., Pfister T., Ricketts T., Robinson L. B., 
Stover R. J., Tucker D., Ward J., Wei M. Z., 1994, SPIE, 
2198, 362 



Wang Z., Chakrabarty D., Roche P., Charles P. A., Kuulk- 
ers E., Shahbaz T., Simpson C, Forbes D. A., Helsdon 
S. F., 2001, ApJ, 563, L61 

Wijnands R., 2006, in Lowry, J. A., Ed., Trends in Pulsar 
Research, Nova Science Publishers, New York, pp. 53-78 

Wijnands R., van der Klis, M., 1998, Nature, 294, 344 

This paper has been typeset from a T^JC/ I^T^^X file prepared 
by the author. 



